Acoustic phonon-limited mobility in two-dimensional M0S2: Deformation potential 
and piezoelectric scattering from first principles 
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We theoretically study the acoustic phonon-limited mobility in two-dimensional n-doped M0S2 
for temperatures T < 100 K and high carrier densities using the Boltzmann equation and first- 
principles calculations of the electron-phonon interaction. The individual deformation potential and 
piezoelectric interactions with the acoustic phonons are extracted using a real-space partitioning 
scheme. Combined with analytical considerations, we find that the piezoelectric interaction in 
two-dimensional (hexagonal) crystals differs qualitatively from the three-dimensional bulk form. 
Owing to the specific form of the piezoelectric interaction in 2D, the two coupling mechanisms are 
found to contribute comparably to carrier scattering throughout the considered temperature range. 
The calculated mobility for screened electron-phonon interaction shows a nontrivial temperature 
dependence both above and below the Bloch-Griineisen temperature. In the low-temperature Bloch- 
Griineisen regime the mobility approaches a fi ~ T -6 (fi ~ T -4 ) behavior for screened (unscreened) 
electron-phonon interaction. Very high intrinsic mobilities in excess of 10 s cm 2 /V s are predicted 
for T < 10 K and high carrier densities (n > 10 11 cm -2 ). 

PACS numbers: 72.10.-d, 72.80.Jc, 81.05.Hd 



I. INTRODUCTION 

Two-dimensional (2D) atomic crystals^ such as 
graphene^ — are promising candidates for future elec- 
tronic applications. Single-layer M0S2 belongs to a 
new family of 2D materials based on transition metal 
dichalcogenides (MX2) which have interesting electronic 
and optical properties^ - — Unlike graphene, 2D M0S2 is 
a regular semiconductor with a band gap of ~1.8 eV£ 
In conjunction with the excellent gate control inherent 
to atomic-thin materials the finite gap makes it a desir- 
able material for various electronic applications. How- 
ever, in spite of the recent progress in sample fabrication 
and transport measurements on gated single to few-layer 
samples ) 1 ' 8 ' 10 : 12- — little is still known about the intrinsic 
carrier properties such as factors limiting the achievable 
mobilities. 

Experimentally, the room-temperature mobility 
of n-type single-layer M0S2 ranges from ~1 to 
^200 cm 2 V -1 s -1 depending on the device struc- 



ture. 
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The higher value obtained for a top-gated 
sample where a high-K gate dielectric was used^ is close 
to the theoretically predicted intrinsic phonon-limited 
value of 410 cm 2 V -1 s -1 for the mobility^ indicating 
that impurity scattering can be strongly suppressed 
by dielectric engineering i±i Other theoretical studies 
have addressed different pertaining issues related to the 
performance of single-layer M0S2 transistors! 18 ' 19 

One important limiting factor for the carrier mo- 
bility of semiconductor-based two-dimensional electron 
gases (2DEGs) is scattering on the acoustic phonons 
of the material. At low-temperatures where scattering 



on optical phonons is suppressed and with weak impu- 
rity scattering, this is among the dominating scatter- 
ing mechanisms In the low-temperature regime, acous- 
tic phonon-dominated transport manifests itself in a 
strong change in the temperature dependence of the car- 
rier mobility once the temperature is lowered below the 
Bloch-Griineisen (BG) temperature Tbg corresponding 
to the thermal energy of an acoustic phonon required 
for full backscattering at the Fermi surface. Explicitly, 
fcjjTjjQ = 2hc\kp, where kp is the Fermi wave vector, 
c\ the sound velocity and &b the Boltzmann constant. 
For heterostructure-based 2DEGs the BG regime is well- 
known ; 21 : 22 and recently transport in the BG regime has 
also been studied in graphene both experimentally^ 3 , and 
theoretically^ - — 

In a 2DEG the Fermi wave vector kp scales with the 
carrier density as y/n and the Bloch-Griineisen temper- 
ature acquires a similar density dependence k^T^o — 
2hc\ \J Ann / g s g v with g s and g v denoting the spin and 
valley degeneracy, respectively.— For single-layer M0S2 
(g v = 2) this results in the following BG temperatures 
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for the transverse (TA) and longitudinal (LA) acoustic 
phonons, respectively, with the carrier density measured 
in units of 10 12 cm -2 . These numbers are on the same 
order of magnitude as those for graphenej— and trans- 
port in the high-mobility BG regime should be achiev- 
able in single-layer M0S2 (and other 2D transition metal 
dichalcogenides). However, the above considerations also 
emphasize the importance of going to high carrier densi- 
ties n > 10 12 cm -2 in order for the BG transition to occur 
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at sufficiently high temperatures where acoustic phonon 
scattering still dominates over impurity scattering^ In- 



deed, extremely large carrier densities n 
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have recently been achieved in 2D samples of graphene 
and M0S2 using advanced electrolytic gatingi 10 ' 23 

In the present work, we study the acoustic phonon- 
limited mobility of n-type 2D M0S2 at low tempera- 
tures (T < 100 K) taking into account both deforma- 
tion potential and piezoelectric scattering. In our pre- 
vious work considering scattering on both acoustic and 
optical phonons, 16 only the deformation potential inter- 
action was taken into account in the coupling to the 
acoustic phonons. We stress, however, that since opti- 
cal phonon scattering dominated the mobility in tem- 
perature range (T > 100 K) considered in Refill the 
conclusions reached there are unaffected by the inclusion 
of piezoelectric scattering. For the temperatures con- 
sidered in this work, scattering on intrinsic intervalley 
acoustic/optical phonons is strongly suppressed and can 
be neglected*^ 

Using a first-principles approach, we here calculate 
both the deformation potential and piezoelectric inter- 
actions in 2D M0S2. Furthermore, an analytic expres- 
sion for the wave-vector dependence of the piezoelectric 
interaction in 2D hexagonal crystals is given. This al- 
lows us to establish the individual coupling strengths of 
the two scattering mechanisms. The calculated transport 
characteristics provide a platform for comparison with fu- 
ture measurements of the carrier mobility in single-layer 
M0S2 which can (i) lead to an experimental verification 
of the theoretical deformation potential and piezoelectric 
couplings reported here^ (ii) reveal to which extent the 
mobility is affected by extrinsic surface acoustic/optical 
phonon o 29 ' 30 which have turned out to be important in 
graphene samples^ - — and (iii) address the importance 
of the interplay between acoustic phonon and impurity 
scattering which results in a complex temperature and 
density dependence of the mobility.— In this context pre- 
vious studies have emphasized the importance of includ- 
ing both the TA and LA phonons in order to obtain good 
agreement with experiment^ 2 . 

The paper is organized as follows. Section |TT] briefly 
summarizes the Boltzmann transport theory. In Sec. IIIII 
the first-principles results for the acoustic electron- 
phonon interaction are presented and the individual de- 
formation potential and piezoelectric interactions are dis- 
cussed. Finally, the results for temperature and density 
dependence of the acoustic phonon-limited mobility are 
presented in Sec. [TV] 



II. BOLTZMANN THEORY 

Two-dimensional M0S2 has a hexagonal lattice struc- 
ture like graphene with the bottom of the conduction 
band residing in the K, if'-points at the corners of the 
Brillouin zone i 35 i 36 The two K, AT'-valleys are perfectly 
isotropic and can to a good approximation be described 



by parabolic bands with an effective mass of m* — 
0.48 m e J£ 

In Boltzmann theory, the drift mobility fj, xx = a xx /ne, 
where a xx is the conductivity, is in the presence of (quasi) 
elastic scattering given by the Drude-like expression^ 7 - 
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where is the energy-dependent relaxation time and the 
energy- weighted average (•) is defined by 



(A) = - /(fe k /9(ek)£k J 4k ( ~^r- 
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Here, n is the two-dimensional carrier density, p = 
g s g v m* /2irh 2 is the density of states in 2D, = 
h 2 k 2 /2m* is the carrier energy, and /(ek) = 
{1 + exp [(^k — Ep)/kBT]}~ 1 is the equilibrium Fermi- 
Dirac distribution function. For a degenerate electron 
gas, only scattering within a thin shell of width k^T 
around the Fermi level is relevant and \i xx ~ eTk F /m* 
applies. 

In the case of acoustic phonon scattering which to a 
good approximation can be treated as a quasielastic scat- 
tering process, the relaxation time for each of the acoustic 
phonons is given b y 16 ' 37 
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where A denotes the acoustic phonon branch index 
(A =TA, LA), #k, k' is the scattering angle and /k = /(ek) 
is understood. The transition matrix element is given by 
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> (5) 



where g£ q is the electron-phonon coupling, e(q,T) is the 
wave vector and temperature dependent static dielectric 
function of the carriers, and hujq\ the acoustic phonon 
energy. The phonons are assumed to be in equilibrium 
and populated according to the Bosc-Einstein distribu- 
tion function Nq\ — N(hiOq\). Under the assumption 
of quasielastic scattering, the phonon energies in the 5- 
functions of Eq. ([5]) are neglected. 

Above the BG temperature where the equipartion ap- 
proximation iV q kBT/hujq 3> 1 applies and with 
screening neglected, the relaxation time for acoustic de- 
formation potential scattering is given by the simple ex- 
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where is the acoustic deformation potential. 38 Here, 
the independence on the carrier energy gives rise to a 
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/i ~ T _1 temperature dependence of the mobility in the 
high-temperature regime. 

The screening of the electron-phonon interaction by 
the carriers themselves is accounted for at the level of 
finite-temperature Thomas- Fermi theory where the static 
dielectric function of the 2DEG is given by 



<q,T) = l 



(7) 
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where qr-p{T) is the finite-temperature Thomas-Fermi 
screening wave vector (see App. [3] for details). 

For the numerical evaluation of the relaxation time in 
Eq. ((4]) quasielastic scattering on the acoustic phonons 
is assumed. It should be stressed that while the phonon 
energies in the (5-functions of Eq. ([5]) are omitted, they 
must be retained in the rapidly varying Fermi functions 
in Eq. (|4j) . This is particularly important in the BG 
regime where the phonon energy becomes comparable to 
the thermal smearing at the Fermi level, i.e. hiOq\ ~ ksT. 



III. INTERACTION WITH ACOUSTIC 
PHONONS 

In the following, the electron-phonon interaction is cal- 
culated from first principles using the density-functional 
based method outlined in RefUa. The interaction with 
the acoustic phonons can be written in the general form 




(8) 



where A is the area of the sample, p is the mass den- 
sity, M£ = (k + q|(5V^A |k) is the matrix element be- 
tween the Bloch states with wave vectors k and k + q, 
and SVq\ is the change in the crystal potential due to 
a phonon with wave vector q and branch index A. The 
couplings in the K, K' valleys are related through time- 
reversal symmetry as | M^ x | = | M_ ' x | . As the hexagonal 
lattice of two-dimensional M0S2 lacks a center of sym- 
metry, charge carriers in single-layer M0S2 interact with 
acoustic phonons through both the deformation potential 
and the piezoelectric interaction. The coupling matrix 
element therefore has contributions from both coupling 
mechanisms, i.e. 



M. 



qA 
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In three-dimensional bulk crystals, the two coupling 
mechanisms are most often assumed to be out of phase, 
i.e. one is real and the other imaginary^ This implies 
that piezoelectric and deformation potential interactions 
do not interfere in lowest-order perturbation theory, i.e. 
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By noticing that the deformation potential is short- 
ranged and the piezoelectric interaction long-ranged 
in nature, we have deviced the real-space partitioning 
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FIG. 1: (Color online) Acoustic electron-phonon couplings 
in the if-valley of the conduction band in two-dimensional 
M0S2. The plots show the absolute value of the first-principles 
coupling matrix elements |M^ pz | at k = K for the TA and 
LA phonons as a function of the two-dimensional phonon 
wave vector q. (Upper row) Deformation potential couplings. 
(Lower row) Piezoelectric couplings. 



scheme outlined in App.|B]to separate out the two mecha- 
nisms from the total electron-phonon coupling. The cou- 
plings to the TA and LA phonons for k = K resulting 
from this partitioning scheme are shown in Fig. Q] with 
the deformation potential and piezoelectric couplings in 
the upper and lower row, respectively^ While the de- 
formation potential couplings have the three-fold rota- 
tional symmetry of the conduction band in the vicinity 
of the K, K' points of the Brillouin zone, the six-fold 
rotational symmetry of the piezoelectric couplings origi- 
nates from the hexagonal crystal lattice. From the first- 
principles calculation of the total matrix elements (not 
shown) we find that the above-mentioned out-of-phase 
property of the two coupling mechanisms only holds for 
the LA phonon. For the TA phonon the two mechanisms 
interfere and must hence be considered together. 

In the following two sections the two coupling mech- 
anisms are analyzed in closer detail. Analytic expres- 
sions for the couplings are introduced and the cou- 
pling strengths are determined from the first-principles 
electron-phonon couplings. 



A. Deformation potential interaction 

The origin of the deformation potential interaction lies 
in the local changes of the crystal potential that accompa- 
nies the atomic displacements due to an acoustic phonon. 
It therefore requires microscopic considerations as the 
one adopted here to determine the interaction strength. 
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In the following we adopt the common approach and 
assume that the deformation potential interaction is 
isotropic and linear in the phonon wave vector, i.e. 

i^ q d ii = s A9 , (io) 

where S> is the acoustic deformation potential.- 8 As this 
is not the case for the true deformation potential cou- 
plings in Fig. [T] which are anisotropic and show a more 
complex g-dependence at shorter wavelengths, the defor- 
mation potential must be regarded as an effective cou- 
pling parameter. In valley-degenerate semiconductors 
where the couplings in the different valleys are related 
through time-reversal symmetry, the effective deforma- 
tion potential of Eq. (|10j) can furthermore be regarded 
as a valley-averaged coupling parameter that implicitly 
accounts for the variation in the coupling between the 
different valleys. We will here use the effective deforma- 
tion potentials for the TA and LA phonons from Refllta 
(see Table [Q. 

It should be noted that the first-principles deformation 
potential couplings contain contributions from umklapp 
processes, i.e. Fourier components of the crystal poten- 
tial involving reciprocal lattice vectors G. As umklapp 
processes are of shorter wavelength (q + G) compared 
to the phonon wave vector q, screening of the electron- 
phonon interaction with e(q) as done in Eq. ([5]) can lead 
to an overscreening of the deformation potential interac- 
tion. However, we are at present not able to quantify the 
relative contributions from normal (G = 0) and umklapp 
processes in further detail. Under the assumption that 
the former dominates the deformation potential interac- 
tion, screening with a long- wavelength dielectric function 
is legitimate. On the other hand, if umklapp processes 
dominate, the individual q + G contributions to the de- 
formation potential interaction should be screened with 
the dielectric function evaluated at q + G. 

B. Piezoelectric interaction 

Piezoelectric coupling to acoustic phonons occurs in 
crystals lacking a center of symmetry and originates from 
the macroscopic polarization that accompanies an ap- 
plied strain . The strength of the interaction is given 
by the piezoelectric tensor e^ here given in Voigt nota- 
tion. 

In App. [C] we provide a derivation of the piezoelec- 
tric interaction in 2D materials based on a continuum 
description. We find that the piezoelectric coupling for a 
2D hexagonal crystal is given by 

\Kx\ = — gxerfc(g<7/2)|^(q)|, (11) 

where en (units of C/m) is the only independent com- 
ponent of the piezoelectric tensor of the 2D hexagonal 
lattice^ eo is the vacuum permeability, erfc is the com- 
plementary error function, a is an effective width of the 



TABLE I: Material parameters for single-layer M0S2 used in 
the present work. Except for the piezoelectric constant which 
has been obtained by fitting the analytic result for the piezo- 
electric interaction in Eq. pip to the first-principles results 
in Fig. [T] all parameters have been adopted from Reffl6l. 



Parameter 

Lattice constant 

Ion mass density 

Effective electron mass 

Transverse sound velocity 

Longitudinal sound velocity 

Acoustic deformation potentials 

TA 

LA 

Piezoelectric constant 
Effective layer thickness 



Symbol Value 

a 3.14 A 

p 3.1 x 10" 7 g/cm 2 

m* 0.48 m e 

cta 4.2 x 10 3 m/s 

c LA 6.7 x 10 3 m/s 

Eta 1.6 eV 

H LA 2.8 eV 

en 3.0 x 10" 13 C/cm 

a 4.41 A 



electronic wave functions, and Ax(q) is an anisotropy 
factor that accounts for the angular dependence of the 
piezoelectric interaction. The latter is given by 

A TA (q) = sin 3 6» q - 3sin6» q cos 2 6l q , (12) 
^LA(q) = cos 3 6» q - 3cos6» q sin 2 q (13) 

for the TA and LA phonons, respectively, — and results 
in the highly anisotropic piezoelectric interaction shown 
in the lower panel of Fig. [TJ The extrema of the piezo- 
electric interaction occur along high-symmetry directions 
of the hexagonal lattice with the directions of minimum 
coupling to the TA phonon and maximum coupling to 
the LA phonon, and vice versa, coinciding. 

The agreement between Eq. (fTTj) and the first- 
principles results for the piezoelectric interaction in Fig.Q] 
is very good (see also Fig. [5] in App. [UJ. Though, for 
q < 2it / L cc \\, where L ce \\ is the size of the supercell used 
for the calculation of the electron-phonon interaction, the 
correct analytic behavior of the first-principles piezoelec- 
tric interaction cannot be captured (see App. [B]). From 
the first-principles results, the piezoelectric constant of 
2D MoS 2 is estimated to be en ~ 3.0 x 10~ 13 C/cm 
(~0.01 e/bohr) which is an order of magnitude smaller 
than the theoretical value for 2D hexagonal BN obtained 
with a Berry's phase approach. 40 Similar Berry's phase 
calculations for 2D M0S2 are required to verify the value 
predicted here. 

The g-dependence of the piezoelectric coupling in 
Eq. (jlip is qualitatively different from the one in 3D 
bulk systems where ~ constant^ In the long- 

wavelength limit where erfc(gcr/2) ~ 1 — qo j\/ 7 2/K, the 
2D piezoelectric interaction in Eq. acquires a lin- 
ear ^-dependence ~ q. Hence, there is no quali- 
tative difference between the deformation potential in- 
teraction in Eq. (|10p and the 2D piezoelectric interac- 
tion in the long- wavelength limit. Under the assumption 
that the linear long-wavelength behavior holds, the high- 
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temperature relaxation time for piezoelectric scattering 
is given by Eq. ^ with the replacement^ 



ene 
eo 



(14) 



where the factor 1/2 = (A\) stems from the angu- 
lar mean of the piezoelectric interaction. This corre- 
sponds to an effective isotropic piezoelectric coupling 
with A x {q) = l/\/2 in Eq. (fTTj) . 

Due to the similarity of the deformation potential and 
piezoelectric interaction in the long- wavelength limit, the 
relative strength of the two interaction mechanisms is 
governed by the ratio (ene/eo)/SA of the prefactors in 
Eqs. (fT0|) and (|TT|) . Since this is on the order of unity 
with the parameters for M0S2 listed in Table HI both cou- 
pling mechanisms must be considered. Furthermore, due 
to the interference between the two mechanisms for the 
TA phonon, the total coupling for the TA phonon differs 
in the K, K 1 valleys for a given direction of the carrier 
wave vector k. The Boltzmann equation must hence be 
solved explicitly in both valleys and the valley-averaged 
relaxation time r = (r K + r A „)/2, where T K/K , denotes 
the total relaxation time in the individual valleys, must 
be used in Eq. ^ for the mobility. 



IV. RESULTS 

In the following we use the Boltzmann equation ap- 
proach outlined in Section [TT] to study the tempera- 
ture (fi ~ T~ 7 ) and density dependence of the acous- 
tic phonon-limited mobility in the temperature regime 
T < 100 K and for high carrier densities 1 x 10 10 - 
3 x 10 13 cm~ 2 . The temperature dependence of the re- 
sistivity p = (ne/i) -1 follows the same power law with a 
change in the sign of the exponent. 

In order to establish the relative strength of deforma- 
tion potential and piezoelectric scattering in 2D M0S2 
we first consider the two scattering mechanisms sepa- 
rately. Figure [5] shows the calculated inverse momentum 
relaxation times due to deformation potential (dashed 
lines) and piezoelectric (full lines) scattering at temper- 
atures T = 10 K and T — 50 K and a carrier density 
of n = 10 13 cm~ 2 . This corresponds to a Fermi level of 
Ep ~ 25 meV and BG temperatures of ~36 K and ^57 K 
for the TA and LA phonons, respectively. The dip in the 
inverse relaxation time at the Fermi level that develops 
with decreasing temperature is a signature of transport in 
the BG regime. In this temperature regime the freezing 
out of short-wavelength acoustic phonons and the sharp- 
ening of the Fermi surface strongly limit the phase space 
available for scattering on acoustic phonons and result 
in a strong suppression of scattering at the Fermi level. 
The figure clearly demonstrates that deformation poten- 
tial and piezoelectric scattering are on the same order 
of magnitude in two-dimensional M0S2. At higher en- 
ergies deformation potential scattering dominates due to 
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FIG. 2: (Color online) Inverse relaxation time for screened 
deformation potential (dashed lines) and piezoelectric (full 
lines) scattering at temperatures T = 10 K and T = 50 K. 
The carrier density is n — 10 13 cm -2 corresponding to a Fermi 
energy of Ef ~ 25 meV and a BG temperature of Tbg ~ 
36 (57) K for the TA (LA) phonon. At T = 10 K the dip 
in the inverse relaxation time at the Fermi energy is a clear 
fingerprint of Bloch-Griineisen physics. 



the different g-dependencies of the two coupling mech- 
anisms. In the relevant energy range around the Fermi 
level, the relative strength of the two scattering mech- 
anisms is found to be weakly dependent on the carrier 
density. 

Due to the highly anisotropic piezoelectric interaction 
and its interference with the deformation potential inter- 
action in the coupling to the TA phonon, the mobility in 
2D M0S2 depends on the direction of the applied elec- 
tric field relative to the orientation of the crystal lattice. 
However, in the following we shall focus on a single ori- 
entation^ 3 - and note that the mobility changes slightly 
with the field oriented along the other high-symmetry 
direction of the hexagonal lattice. Along the chosen di- 
rection of the applied field, we find that the interference 
in the coupling to the TA phonon cancels and the de- 
formation potential and piezoelectric interactions can be 
treated separately. 

Before we consider the mobility in the presence of car- 
rier screening, we briefly summarize the temperature de- 
pendence of the mobility in the absence of screening. 
This is relevant at low carrier densities and for deforma- 
tion potential couplings dominated by umklapp processes 
where screening by the carriers is less efficient and to a 
good approximation can be neglected. Above the BG 
temperature, the relaxation time for deformation poten- 
tial scattering in Eq. ^ gives rise to a linear (7 = 1) be- 
havior while piezoelectric scattering produces a slightly 
sublincar (7 < 1) behavior. The total mobility follows 
the latter. At low temperatures T <§C Tbg, we find that 
the temperature dependence converges towards a 7 = 4 
behavior in all cases due to the identical long-wavelength 
limits of the deformation potential and piezoelectric cou- 
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FIG. 3: (Color online) Acoustic phonon-limited mobility vs 
temperature for screened electron-phonon interaction. (Up- 
per) Mobility vs temperature for different carrier densities. 
(Lower) Temperature dependence of the exponent 7 in fi ~ 
T -7 for the same set of carrier densities. 



plings in Eqs. (JTOJ) and (|TT|) . 

Figure [3] summarizes the temperature dependence of 
the calculated acoustic phonon-limited mobility for car- 
rier densities 1 x 10 10 - 3 x 10 13 cm -2 corresponding 
to BG temperatures up to -62 K (-99 K) for the TA 
(LA) phonon. Both the mobility (upper) and the ex- 
ponent 7 = — d log n/d log T of its temperature depen- 
dence (lower) are shown. In general, the interplay be- 
tween BG physics and the quantum-classical crossover 
from a degenerate to a nondegenerate carrier distribu- 
tion at T — Tp w 29tj K (with n measured in units of 
10 12 cm -2 ) results in a nontrivial temperature and den- 
sity dependence of the mobility. At the lowest densities 
the characteristic temperatures Tbg and Tp are compa- 
rable while Tp > Tbg for n > 10 12 cm -2 . The crossover 
to the BG regime at T — Tbg is therefore more clear for 
the higher carrier densities. 

At temperatures T > Tp above the quantum-classical 
crossover, the mobility shows a 1 < 7 < 4 dependence for 
all densities which originates from the temperature de- 
pendent the dielectric function. The characteristic non- 
monotic behavior of the exponent in this regime shifts 
to higher temperatures with increasing density. A more 
correct description of the dielectric function might give 
rise to minor quantitative modifications of this tempera- 



ture dependence (see e.g. Refl34l). The unscreened high- 
temperature behavior /i — T -1 of the mobility is only ap- 
proached for the lowest carrier density in the considered 
temperature range. The expression = l/^o + oT 
often used to fit experimental mobilities^ where /xo is 
the residual mobility due to e.g. impurity scattering and 
a is the density-dependent slope of the temperature de- 
pendence, is therefore not expected to apply. 

In the low-temperature BG regime T < Tbg, a 
stronger 7 > 4 temperature behavior appears. The non- 
monotonic behavior of the exponent in the BG regime is 
found to be present also in the absence of screening and 
can therefore not be attributed to screening effects. As 
expected from Eq. (JT|), the characteristic features in the 
temperature dependence of the exponent also here shift 
to higher temperatures with increasing density. We find 
that the mobility approaches a 7 = 6 limiting behav- 
ior at T <C Tbg- The same holds for the mobilities due 
to separate deformation potential and piezoelectric scat- 
tering. This is in agreement with the analytic considera- 
tions for deformation potential scattering on 2D phonons 
of Ref 34. Because of the specific form of the piezoelec- 
tric interaction derived here and the 2D nature of the 
phonons, this differs from the usual low-temperature de- 
pendence for a 2DEG with scattering on bulk 3D phonons 
in heterostructures where the temperature dependence is 
predicted to saturate at 7 = 5 and 7 = 7 for screened 
piezoelectric and deformation potential interaction, re- 
spectively jSii^ As the calculated low-temperature lim- 
its of the mobility are only realized deep inside the BG 
regime T <C Tbg, they m ay, however, be of limited prac- 
tical relevance. 

The low-temperature mobilities at T < 8 K show a 
general increase as a function of the carrier density. At 
T > 8 K the change in the screening properties follow- 
ing the quantum-classical crossover makes the density de- 
pendence more complicated. The density dependence of 
the mobility is shown explicitly in Fig. 2] for T = 4 K 
and T = 20 K. At densities corresponding to a Fermi 
temperature Tp > T slightly larger than the tempera- 
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FIG. 4: (Color online) Acoustic phonon-limited mobility vs 
carrier density for screened electron-phonon interaction and 
temperatures T = 4 K and T = 20 K. 
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ture, the increasing low-density behavior is replaced by 
a plateau with a slowly decreasing density dependence. 
Only for degenerate carrier distributions n > 10 12 cm" 2 
at T = 4 K is this followed by an increasing density de- 
pendence. The increase in the mobility with the carrier 
density in the degenerate low-temperature regime can be 
explained by the corresponding increase in the Fermi ve- 
locity. Other factors affecting the mobility such as the 
level of screening and phonon scattering do not vary 
significantly with the density in this regime. Together 
with the general high level of screening in 2D M0S2 (see 
App. [A"]l . this results in very high intrinsic mobilities in 
excess of 10 s cm 2 V -1 s _1 at T < 10 K and high carrier 
densities. 



V. CONCLUSIONS 

In this work we have combined analytic and first- 
principles calculations of the electron-phonon interaction 
with a numerical solution of the Boltzmann equation to 
study the acoustic phonon- limited mobility in 2D M0S2. 
Our findings are also relevant for the carrier mobility in 
other 2D piezoelectric materials. The key results are, (i) 
the analytic expression in Eq. for the piezoelectric 
interaction in 2D hexagonal materials which is qualita- 
tively different from the one in 3D bulk system and which 
results in similar temperature dependencies of deforma- 
tion potential and piezoelectric scattering, (ii) a method 
to separate out short-range and long-range contributions 
to the electron-phonon interaction in first-principles cal- 
culations, (iii) fi ~ T-T with 7 - 6 (4) at T < T BG for 
screened (unscreened) electron-phonon interaction, and 
(iv) an intrinsic mobility that despite the large effective 
mass (m* = 0.48) of the conduction band in 2D M0S2 
increases beyond 10 s cm 2 V -1 s _1 at T < 10 K and high 
carrier densities. 
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Appendix A: 2D Thomas-Fermi screening 



where 



qTF{T) = g T F 



1 — exp 



-E F 



(A2) 



is the temperature dependent screening wave vector, 
<Ztf = e 2 /(2eo)/Oo the zero-temperature TF wave vector, 
Po = .9sfft)W*/27r/i 2 the constant density of states in 2D, 
and Ep = n/po the Fermi level. This reproduces the 
correct RPA and Debye-Hiickel results 



Qtf(T -> 0) = q TF , 



Qtf (T -> 00) 



ne 



2e fc B T 



(A3) 
(A4) 



for the screening wave vector in the low and high- 
temperature limits, respectively. 

Thomas-Fermi screening suffers from the shortcom- 
ing that it overestimates the screening properties of the 
2DEG at q > 2kp where there is no phase space available 
for screening processes. RPA corrections to the dielec- 
tric function are required to cure this problem.— How- 
ever, considering quasi-elastic scattering processes where 
q < 2kp, TF theory provides a reasonable approximation 
to the dielectric function. 

Due to the density independence of (?tf, zero- 
temperature screening in 2D decreases with increasing 
density (Fermi wave vector kp) when scattering on the 
full Fermi surface is possible. The screening efficiency is 
controlled by the dimensionless ratio q s = <7tf/^f which 
in 2D is given by 



(A5) 



For q s 3> 1 all scattering processes at the Fermi surface 
(q < 2&f) are strongly screened. In the valley-degenerate 
conduction band of 2D M0S2, q s ~ 2 x 10 2 /y / n with n 
measured in units of 10 12 cm -2 . The resulting high level 
of screening is inherent to valley-degenerate 2D semicon- 
ductors with large effective masses. 

From the above considerations we conclude that 
screening in 2D M0S2 is very pronounced and reduces 
the inverse relaxation time in Eq. §Q by several orders of 
magnitude for the carrier densities considered here. How- 
ever, it should be noted that the effect could here be over- 
estimated due to the neglect of the form factor from the 
^-dependence of the electronic wave functions^ First- 
principles calculations of the RPA response function (see 
e.g. Refl45l) for 2D M0S2 could help to shed further light 
on this issue. This is, however, not expected to change 
the level of screening by orders of magnitude. 



In the Thomas-Fermi (TF) approach to screening, the 
finite-temperature dielectric function of a 2DEG is in the 
long- wavelength limit given by 



e(q,T) = l + 



(Al) 



Appendix B: Real-space partitioning scheme for 
separation of the short and long-range part of the 
electron-phonon interaction 

In first-principles calculations of the electron-phonon 
interaction both short-range (deformation potential) and 
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FIG. 5: (Color online) Schematic illustration of the LCAO 
supercell matrix involved in the calculation of the coupling 
Mk q in Eq. (|B1|) . The square lattice indicates the unit cells 
of the crystal lattice. The real-space cutoff r cut measured 
from the position of the atomic site where the gradient of 
the potential is taken, is used to separate out a short-range 
and a long-range contribution to the coupling matrix element. 
Matrix elements involving LCAO orbitals located beyond the 
cutoff as the one shown are defined as long-ranged. 

long-range (piezoelectric and Frohlich interaction) are in 
general included in the coupling given in Eq. ©. How- 
ever, due to their different origin it may sometimes be 
desirable to consider them separately. In the following 
we provide a real-space partitioning scheme to separate 
out the short-range and long-range contributions to the 
coupling matrix element M£ . 

In the present work the electron-phonon interaction 
has been calculated with an LCAO based supercell 
method. 16 Within this approach the coupling matrix el- 
ement due to the phonon-induced change SVq\ in the 
crystal potential is given by 

M k A q = (k + q|JV qA (r)|k) 

= ^E< c ,E elk - (R "" R '" K4q - Rm 

ij mn 

x viRjeqA • VoV(r)|jR„), (Bl) 

where the last equality follows by expanding the Bloch 
states and phonon-induced change of the potential in an 
LCAO basis and atomic displacements, respectively, and 
subsequently exploiting the periodicity of crystal lattice. 
Here, |iR m ) denotes an atomic orbital on atom i of the 
primitive unit cell displaced from a reference cell at Ro 
by the lattice vector R TO , Ci is the LCAO expansion coef- 
ficients, e qA is the mass-scaled phonon polarization vec- 
tor, and VoV is the gradient of the crystal potential with 
respect to atomic displacements in the reference cell. 

The central quantity of the partitioning scheme is the 
LCAO supercell matrix of the atomic gradients in the 
last line of Eq. (|B1[) . The real-space structure of its ma- 
trix elements is illustrated schematically in Fig. [5j The 
partitioning scheme consists in splitting up the summa- 
tion in Eq. (|B1I) into two contributions, (i) a short-ranged 
part which neglects all matrix elements involving LCAO 
orbital beyond a chosen real-space cutoff r cut , and (ii) 
a long-ranged part which includes the remaining matrix 
elements. In general, the real-space cutoff r cut must be 
chosen small enough that the short-range part does not 



include contributions from truly long-range effects in the 
relevant range of phonon wave vectors, i.e. 7r/r cut > q max 
where Qmax is the maximum phonon wave vector of inter- 
est. At the same time, the cutoff cannot be chosen too 
small that short-range effects are cut off. The finite real- 
space range of the electron-phonon interaction inherent 
to a supercell approach naturally sets a lower limit q m i n 
for the magnitude of the phonon wave vector q at which 
the long-ranged interactions be obtained reliably. It is 
given by = 27r/L ce n where L ce n is the size of the 
supercell (measured as the diameter of a sphere that can 
be contained within the supercell). 

The real-space partitioning scheme outlined here can 
be applied in other first-principles calculations of the 
electron-phonon interaction based on e.g. Wannier func- 
tions 4 s - 



Appendix C: Piezoelectric interaction in 2D 
hexagonal crystals 

Piezoelectric interaction with acoustic phonons ap- 
pears in crystals which lack a center of symmetry. In this 
case, the displacement field u q (r) of the acoustic phonons 
leads to a polarization P of the lattice given by2£ 

/' X/'- 4 '"- hj' k = x ^y ( C1 ) 

jk 

where 

*-K£ + &) (C2) 

is the strain tensor and e is the tensor of piezoelectric 
moduli having symmetry &tjk = e i,kj- I n 2D materi- 
als the piezoelectric coupling has units of C/m (C/m 2 
in 3D) — the displacement field can be thought to have a 
normalized ^-profile with units of m -1 . For a 2D hexag- 
onal lattice with a basis there is only one independent 
piezoelectric component en (Voigt notation) which is re- 
lated to the other nonzero components as^S 

en = -ei2 = -e 26 . (C3) 

Here, the lattice vectors of the hexagonal lattice has been 
chosen as a^a = a(v3/2, ±1/2) where a is the lattice con- 
stant. For two-dimensional h-BN en ~ 1CP 12 C/cm4°- 

In the following we derive an analytic expression for 
the piezoelectric interaction in 2D materials. Consider- 
ing an isolated 2D material sheet, the in- plane atomic 
displacements due to an acoustic phonons are in a con- 
tinuum description represented by the 2D displacement 
field 

u qA (r) =e A e^ r / q (z)Q qA (C4) 

where q = (q Xl q y ) is the two-dimensional phonon wave 
vector, e A is a unit vector describing the polarization of 
the phonon A (to a good approximation independent on 
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FIG. 6: (Color online) Piezoelectric coupling in a 2D hexago- 
nal crystal. The plots show the absolute value of the coupling 
matrix element Mp Z (q) in Eg. ICTol for the TA (left) and LA 
(right) phonons as a function of the phonon wave vector q. 
The parameters for M0S2 in Tab. |T] have been used. 

q in the long- wavelength limit) ^ / q is the z-profile of 
the displacement field in the direction perpendicular to 
the material sheet, and Qq\ = ^h/2Apuj cl \(al lX + a-q\) 
is the vibrational normal coordinate. Since we are here 
interested in the coupling matrix element M^ q in Eq. (|8]) 
the normal coordinate will be left out in the following. 

The components of strain tensor follow directly from 
Eq. (|C2]) as 

(■xx = iqxe\, x e l ' 1 ' r fq(z) (C5) 
e X y = £yx = - (q y ex,x + q x e\,y) e iq ' r / q (z), (C6) 
and the associated polarization in Eq. (|C1[) is given by 

P q,xi r ) = e H i e xx ~ tyy) 

= ie n (q x e\. x - q v e\, y ) e lq ' r / q (z) (C7) 
O r ) = -2e lie ^ 

= -ien (q y e\,x + qxe\,y) e iq ' r / q (z). (C8) 

The potential <f> resulting from the piezoelectric polariza- 
tion field is given by Poisson's equation — eoV 2 (/>(r) = p 
where p = — V • P is the polarization charge. Since we are 
considering an isolated material sheet, the only bound- 
ary condition that applies is 4> — > for z — > ±00. Fourier 
transforming in all three directions, we find 

eo(q 2 + k 2 )^(k) = p q / q (fc), (C9) 

where k is the Fourier variable in the direction perpendic- 
ular to the plane of the layer and the Fourier components 
of the polarization charge are given by 

PqA = -iq • P q 

= e ii (q x e\,x - q^Kx - 2qxq y e\, y ) 

= e liq 2 A x (q) (CIO) 



where the angular dependencies have been collected in 
the anisotropy factors A\. It is given by 

At a (q) = sin 3 q - 3 sin6» q cos 2 q (Cll) 

A LA (q) = cos 3 (9 q -3cos(9 q sin 2 (9 q , (C12) 

for the TA (Sta -L q) and LA (Sta || q) phonons, respec- 
tively. 

The z-dependence of the potential is given be the in- 
verse Fourier transform which yields 

r 2 

<M*) = — Mk) / dk e lkz Uk)^— 
eo J q + « 

= £ii^ A (q)ge-«l-l, (C13) 
(a 

where in the last equality, a (5-function z-profile with 
fq(k) = 1 has been assumed. For atomically thin ma- 
terials, this should be a reasonable approximation. 

The coupling matrix element in Eq. (|5J| follows from 
the overlap with the envelope x( z ) of the electronic Bloch 
functions, 

M*M = jdz x*(*)<M*)Xv2)- (C14) 

Assuming, for simplicity, a Gaussian envelope the long- 
wavelength limit of the coupling is given byi2. 

M*M = e -^q x erfc(«<r/2)A A ($. (C15) 

where a is the effective width of the electronic envelope 
function. Contrary to the 3D case where the piezoelec- 
tric coupling is independent on the coupling here 
depends on the magnitude of the phonon wave vector. 
The absolute value of the coupling matrix elements are 
shown in Fig. [5] The six-fold rotational symmetry intro- 
duced by the anisotropy factors in Eq. (IClip stems from 
the hexagonal lattice. 

A more detailed treatment of the electronic envelope 
function is not expected to change the g-dependence of 
the interaction in Eq. (IC15I) qualitative. Indeed, the good 
agreement with our first-principles calculations shown in 
Fig. [U confirms this expectation. The introduction of dif- 
ferent boundary conditions for the electrostatic potential 
relevant for e.g. substrate-supported materials sheets is 
also not expected to change things qualitatively. In any 
case, the strength of the piezoelectric interaction is given 
by the piezoelectric constant en of the material and is 
not affected by such details. 
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